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ABSTRACT 

We present Herschel images from the HOBYS key program of the Eagle Nebula (Ml 6) in the far-infrared and sub-millimetre, using the PACS 
and SPIRE cameras at 70 jim, 160 jim, 250 jim, 350 jim, 500 jim. Ml 6, home to the Pillars of Creation, is largely under the influence of the nearby 
NGC 661 1 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved 
by the OB cluster. We investigate the heating effect of NGC 6611 on the entire Ml 6 star-forming complex seen by Herschel including the diffuse 
cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of 
Ml 6 with respect to the nebula, its surrounding environment, and the NGC 6611 cavity. The dust temperature and column density maps reveal a 
prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC 6611 cluster, as well as a 
northern filament which extends around and away from the cluster across a potential ridge of star formation. The dust temperature in each of these 
filaments decreases with increasing distance from the NGC 6611 cluster, indicating a heating penetration depth of ~ lOpc in each direction in 3 - 
6x10^^ cm~^ column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming 
sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. 
Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which 
was destroyed by the creation of the OB cluster and its subsequent winds and radiation. 

ISM: filaments - ISM: structure - ISM: dust, extinction - Stars: early-type - Stars: pro- 



Key words. ISM: individual objects (M16, NGC 6611) - 
tostars 



1. Introduction 

Forming in dynamic and ever-changing environments, high- 
mass stars are complex entities which, through radiation and 
stellar winds, heavily impact their natal environments, possibly 
even triggering the formation of new stars and planetary systems. 
A general rarity of candidates and their rapid evolutions mean 
that it is difficult to characterise the early e volutionary stages of 
high-mass stars (cf. lGaray & Lizanoll 19991) . 

The He r schel Space Observatory (launched 2009; 
iPilbratt et al.l l2Q10l) has already had a significant impact 
on star formation studies, owing primarily to the large key 
programs devoted to this field, specifically the Herschel iniaging 
study of OB Young Ste llar Objects (HOB YS ;[ Motte et al.l l2'0ia) 
and Gould Beh Survey (Andre et all201 0). as well as the Hi-Gal 
survey (Molina ri et al. 2010 ). These programs target the natal 
environments of bu rgeoning young stars, such as interstella r 

filaments and ridges (lArzoumanian et al ] i2QTTtiHnr^rd]i2QTTh . 

as well as source s embedded within them, e.g., prestellar cores 
and protosta rs (Bontemps et al.l |2QTqI iHennemann et al. 1 120101: 
iGiannini eTal. 2012), and investigate the impact of OB stars o n 
cloud formation (IZavagno et al.l |20TqI: ISchneider et al.l l2010l) . 



The primary HOBYSQ driver is to compile an unbiased census 
of high-mass star progenitors, and their natal ridges, in th e 
nearest (< 3 kpc) star-forming complexes (iMotte et al.l l2010l) . 
one of which is Ml 6. 

The Eagle N ebula (Ml 6) is a young (l-3xl0^yr; 
iHillenbrandetaD 119931) . active high-mass star-forming region 
in the constellation of Serpens ~ 2kpc from the Sun. Ml 6 was 
first observed by de Cheseaux in 1745 and then rediscovered 
by Messier in 1764 who was the first to comment on its neb- 
ulosity. The Europ ean Space Agency's (ESAs) Infrared Space 
Observatory (ISO: ICesarsky et al]|l996h revealed the presence 
of unusually hot dust in the core of the nebula dPilbratt et al.l 
l1998, T colour ~200K). But it w as the Hubble Space Telescope 
images of its ionised gas ( Hester et al.lll996l) . better known as 
the 'Pillars of Creation', that has made Ml 6 iconic and one of 
the most observed objects in the sky. 

Ml 6, and especially these pillars, has since been observed 
at many wavelengths spanning the electromagnetic sp ectrum, 
as well as in many ded i cated spectral line studies (e.g.,|Poun( 



1998; White et al 



y aea i catea spectr al lin e stuaies (e.g.,irounc 
I1999I: lAllen etanfT999l: lUrquhart et al.ll2QQ3l 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

Figs. IB. II & ICll are only available in electronic form at 
khttp : //www. aanda. org 



Linsky et al.il2007l to name a few). More recently, M16 was ob- 
served in the mid-infrared with IRAC (3.6-8.0)Um) as part o f 
the Spitzer GLIMPSE legacy program (llndebetouw et aLlboOTl) . 
and in the far-infrared with MI PS (24 iim & 70 iim ) as part of 
the Spitzer MIPSGAL program ( Elagev et all201 li) . The former 



http ://w w w. herschel . fr/cea/hoby s/en/ 
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15 pc 





Fig. 1. Three-colour image of M16 using 70yum(blue), 160yum(green), 250yum(red). Left: Image of the entire M16 complex with 
well-known regions from the literature, as well as regions identified (the arch and ridge features) and discussed in this paper, 
indicated on the map. Right: Zoom to the central region of Ml 6 (see left image). 



survey targeted the warm dust toward this region and revealed 
its (large) population of young stellar objects (YSOs) that these 
authors deemed protostellar candidates. The latter paper exam- 
ined the dust emission and composition, such as the grain type 
and size, in this region, and concluded that a supernova remnant 
may be responsible for heating the inner portion of Ml 6. 

Adjacent to the Eagle Nebula, and responsible for 
heating and ionising the Nebula, is the young open 
star cluster NG C6611. At a distance of ~ 1.75- 2 kpc 
(iGvaramadze & Bomans 2008), NGC 6611 contains four early- 
type O stars (~ 2 - 3 Myr; Duft on et al .112 006), an d is also associ- 
ated w ith the cluster of a B2.5I star (~ 6 Myr; ^Hillenbrand et al.l 
Il993h . The t otal stellar mass of th e NGC 6611 cluster is 
~2xl O^Mo ('Wolff' etal." '2007"). See Gvaramadze & BomansI 
(|2008|) and Belikov et al. ( 2000) for more information on the 
NGC 6611 cluster. 

We present here unprecedented high-resolution high- 
sensitivity Herschel observations, spanning the far-infrared and 
submillimetre regime, of the Eagle Nebula. These Herschel data 
provide unique access to this far-infrared regime at much greater 
spatial resolution, and additionally cover a greater area around 
the central NGC 6611 cluster and the Pillars of Creation, than 
previous (space) observations and studies. Herschel allows us to 
probe the cold component of the dust, complementary to Spitzer, 
and to deepen the study of warm dust (through the 70 jim band) 
in this region. Here, we present a structural analysis of the en- 
tire Ml 6 molecular cloud complex, with respect to the well- 
known centra l OB cluster. A companion paper by IWhite et aP 
(12012, in pre p) examines the protostellar content of the Ml 6 re- 
gion as detected by Herschel. 

2. Observations and data reduction 

Ml 6 was observed on 2010 March 24 and 2 010 September 11 - 
12, as part of the HOBYS key program dMotte et al.l l2010h . 



In March, the SPIRE (250//m, 350//m, 500 fim: [Griffin et al.l 
2010) camera was inadvert ently turned off and thus only the 
shorter PACS (70, 160yum: iPoglitsch et al.ll2010l) wavelengths 
were mapped during this period. In September, the parallel- scan 
mode of Herschel was used to map simultaneously with both in- 
struments at the aforementioned five bands, using the slow scan 
speed (20'7s). For these September observations, the SPIRE 
camera was aligned with the PACS data taken from March, al- 
lowing 100% common area coverage (normally, the positional 
offset between the PACS and SPIRE instruments results in a 
small (200 spatial offset between the PACS and SPIRE maps), 
as well as an extended PACS coverage to the east and south (see 
Fig. lA.li top). A total area of ~ 1.5 deg by 1.5 deg was mapped 
using two orthogonal scan directions. 

Each of the March and September data were reduced 
with version 5.0.1975 of the H erschel Interactive Processing 
Environment (HIP^; Ottll2010l) adopting standard steps of the 
default pipeline to Level 1 products including calibration and 
deglitching. To improve the baseline subtraction, data taken dur- 
ing the turn-around of the telescope were used. Calibration of the 
PACS data has been found to be within 5%, whilst SPIRE cali- 
bration is within 10% for all bands (see the respective observers' 
manuals). 

Maps were produced using the HIPE Level 1 data and v7 
of the Scanamorphos software packag e^ which perfo rms base- 
line and drift removal before regriding (lRousselJl2012l) . The two 
sets of PACS data, from March and September, were combined 
within Scanamorphos using the '/galactic' option, to produce a 
single map of the entire PACS coverage. For simplicity, we here- 
after use PACS data to refer to the combination of the two indi- 
vidual PACS observations, rather than referring to one particular 
data set. The individual Herschel images at each wavelength can 



^ http : //herschel . esac . esa. int/HIPE_download. shtml 
^ http://www2.iap.fr/users/roussel/herschel/index.html 
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be found in Fig. lA.U including the highest resolution 70 fim im- 
age which traces the hot dust and protostars in this region. The 
three colour image of the the entire Ml 6 complex, as well as a 
zoom to the central region around the OB cluster is given in Fig. 

m 



3. Structural analysis 



The Herschel maps of the Eagle Nebula show bright emission 
throughout, particularly around the NGC 661 1 cluster (see Fig.[T] 
and Eig. lA.lL This bright emission covers the Pillars of Creation 
and the SpireQ east of the cluster and a possible ridge (see Section 
13.31) of material to the north-west. 

Immediately to the north of this ridge is a prominent arch- 
like structure (see Fig.[T] left). While weak diffuse emission can 
be seen toward this arch-like structure in the infrared (8 yum; 
see Eig. 5 of llndebetouw et al. 2007), it is much more visu- 
ally striking and associated with a clear density enhancement 
in our Herschel images e.g., Eig. [TJ right (see also Eig. [2] and 
Eig.[3]below). To the s outh-west of the ridge is th e massive YSO 
(MYSO) identified by llndebetouw et al.l (I2007L also known as 
IRAS 18152-1346) as the most luminous YSO (1 x 10^ L©) in 
the region. The Herschel images of Ml 6 indicate many other 
sites of star formation, such as small clusters of sources in the 
south-east of the map, which are more easily identified in the 
higher resolution PACS data (see Eig. lA.ll top). 

The three-colour composite image of this region (Eig. 
[TJ shows a very distinct temperature gradient away from 
the NGC 6611 cluster. This temperature gradient is reminis- 
cent of that seen in ot her HOBYS regions such as RCW 
120 (IZavagno et al1l2010i) and the bubble Hii region study of 
I Anderson et al.l (l2012l)" 



^The Spir e is a w ell known region in the Eagle Nebula. 
|http : //en . wiki pedia.org/w iki/Mess ier_ 16, 



See, 



Fig. 3. Three colour composite of the column density (red), dust 
temperature (blue) and 70 yum (green) images. The white circle 
indicates the location of the NGC 661 1 cluster. The extent of the 
influence of the photo-dissociation region, created by the pres- 
ence of the OB cluster, is clearly highlighted in green at the cen- 
tre of the map. Limb-brightening efl'ects occur here around col- 
umn density peaks and on many low column density filaments 
around the OB cluster. 



3.1. Column density & dust temperature maps 

The dust temperature and column density N//2) maps of Ml 6 
were drawn by fitting, pixel-by-pixel, spectral energy distribu- 
tions (SEDs) using a modified blackbody model (iHiU et al.l2009l 
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l20ld) . following the procedure outlined by iHill et al.l (l2Qllh . 
These data were convolved to the resolution of the 500 //m band 
(37^0, and the zero offsets, which were determined from compar- 
ison with Planck and IRAS following the procedure detailed by 
iBemard et al.l (l2010h . were applied to the individual maps prior 
to fitting. Only the four longest Herschel wavebands, which trace 
cold dense material, were used to make the column density and 
dust temperature rnaps (Fig, p]), a r nethod consistent within th e 
HOBYS program (iHill et al.ll201 ll: iNguven Luong et aLHSoTl) . 
The 70 yum emission is excluded from the fits as it is not tracing 
the cold dust that we are most interested in, and is likely trac- 
ing small grains in hot PDRs instead (cf. Fig. O. Furthermore, 
excluding the 70 jim emission ensures that the SED fit is more 
robust, since the lOfim flux does not completely fit a sin- 
gle temperature greybody model dominated by longer wave- 
length fluxes. We used a dust opacity law in which the dust 
opacity per unit mass column density (/Cy), is given by Ky = 
0.1 (y/1000 GUzf cmVg (e.g. iMotte etal][2Q10) . and a dust 
spectral index of J3 = 2. The lat ter was fixed in line with the lit- 
erature (iGoldsmith et al.l [79971 and references therein) and be- 
cause we do not have longer wavelength (sub)millimetre data 
to constrain JS accurately. The quality of a SED fit was assessed 
using x'^ minimisation. Most of the individual single tempera- 
ture SED fits comprising the dust temperature and column den- 
sity maps peaked around, or longward of, the 160yum flux - the 
shortest wavelength data used in the fit. The observed rms of cir- 
rus noise in the derived column density map is -2x10^^ cm"^ 
(Ay ~ 2 mag). Note however, that the column density and dust 
temperature are line of sight measures. 

The median temperature throughout Ml 6 is ~18K, with 
slightly cooler temperatures in the northern part of the Herschel 
map, which is furthest from the Galactic Plane. The eastern por- 
tion of M16 is dominated by a dense filament (3 - 7 x 10^^ cm"^) 
situated directly east of, and running parallel to the Pillars of 
Creation. This 'Eastern Filament' appears as a cool temperature 
depression (~16 K; Fig.O right, see Section |4Jl) with respect to 
the warmer surrounding material. 

Figure [3l a three colour composite of the column density, 
dust temperature and 10 /dm image, highlights the region af- 
fected by strong PDR emission, i.e., the interface between the 
cloud and the cavity swept out by the NGC 661 1 cluster. The in- 
clusion of the 70yum emission here highlights clear examples 
of limb brightening, where the filamentary structures around 
the cavity have heated edges at the ionisation front created by 
the NGC 6611 cluster. Furthermore, this Figure demonstrates 
the strong correspondence between high column density and 
low temperature, especially with respect to the Eastern Filament 
(clearly visible here in red). 

3.2. Filamentary structure 

Recent Herschel observations have revealed the ubiquity 
of interstell a r filaments, especially in star-forming regions. 
lAndre et al.l (l2010l l201ll) proposed that in low-mass star- 
forming regions only those filaments above an Ay of ~ 7 mag 
are gravitationally unstable (supercritical) and thus capable of 
forming prestellar cores in their interiors. In fact, a large por- 
tion of the molecular complexes observed as part of HOBYS 
are above this star-forming threshold (e.g.. [Hill et al.l 1201 ll: 
iNguven Luong et al.|[20TllSchneider et al.ll20m 

To take a census of the filaments in the M16 star-forming 
comp lex, the DisPerSE algorithm (Sousbie 2011; Sousbie et al. 
I2OIII) was applied to the column density map. DisPerSE works 
with Morse theory and the concept of persistence to identify 




18^20"'O0' 00' 
Right Ascension (2000) 



Fig. 4. Filaments (cyan) in M16 detected by DisPerSE overlaid 
on the column density map. The Eastern Filament is shown here 
in green, whilst the Northern Filament is in magenta. The po- 
sition of the NGC 6611 cluster is depicted by the yellow star. 
The temperature minimum of the Eastern Filament and the end 
of the temperature gradient for the Northern Filament, seen in 
Fig. [21 are indicated by lines perpendicular to the respective fil- 
aments. The possible Western Filament (see Section 15.11) runs 
parallel to the Eastern Filament, and the possible continuation to 
the Northern Filament is indicated by a dashed line. 

topological structures (see lSousbiell201ll for a detailed explana- 
tion of these concepts) such as filaments and ridges, and connect 
their saddle-points with maxima by integral lines. This algorithm 
has already been successfully applied by our group t o trace fil- 
amen t ary structures in sta r-forming regions (e.g., see tiill et al.l 
1201 ll: iPeretto etal1l2012h . The crest points of the filaments, as 
identified by DisPerSE are presented in Fig. |4] 

The southern half of Ml 6 is dominated by two main fila- 
ments which run roughly parallel to each other near the centre 
of the Herschel map, before converging to a solitary filament 
near the southern edge of the map. The eastern one ('Eastern 
Filament, see Section |4J1 Fig.|4]) is much more prominent (by a 
factor of > 2 with respect to column density, see also Fig.|2l left) 
than the western one. The Pillars of Creation lie on a filamenljj, 
which joins the aforementioned Western Filament (see Fig.|4]). 

The northern half of the Ml 6 region, just north of the central 
NGC 661 1 cluster, displays a complex and clustered network of 
filaments. The most coherent filament in this region starts close 
to the NGC 6611 cluster, connects to the ridge and thereafter 
heads north-west (see Fig.|4]). We hereafter call this filament the 
'Northern Filament' - see Section l42l 

Interestingly the clear 70yum filaments, seen to the south- 
east (see Fig. lA 11 top-left), which highlight PDR or low column 
density filaments, are generally perpendicular to the main SE- 
NW direction of the cloud and the major filaments (see Fig. lA. 11 

^ Note that DisPerSE identifies a single filament for thee pillars, from 
the column density crest points. This is simply a limitation of the col- 
umn density map, as at 37'' resolution only the densest (western most) 
pillar is detected (see Fig-E) right). 
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top-left). These 70yum emission features are not associated with 
point like emission of YSOs, but rather to heating/UV ionisa- 
tion of cloud edges, such as limb brightening/PDRs. These fila- 
ments are similar to the striatio ns seen in regions such as Taurus 
(iPalmeirim et al.ll2012. in prepi) . 

3.3. The arch and a potential ridge 

In the individual Herschel images of the Eagle Nebula there 
is clear evidence of a arch-like structure to the north-west of 
the NGC6611 cluster (Eig. □ and Eig. This arch-like 

structure is actually composed of two arches, an inner warmer 
(~ 17-20K) arch, and an outer high column density (-2.0- 
5.8x10^^ cm"^) arch. These two arches may in fact be tracing 
different parts of the same shell-like structure or bubble, as dis- 
cussed below in Section [5Jl These structures are more promi- 
nent at the shorter PACS wavelengths. Due to its tendency of 
tracing density crests, DisPerSE only traces the column density 
arch (see, Eig.©. 

Since the arch app ears weak at 8yum (with IRAC; see 
llndebetouw et al.l [20071 Eig. 5), and not as clearly as in our 
Herschel images (see for example Eigs. [Tl and lA. 1 1 and Section 
O, it is likely composed of mostly cool mat erial. At the base of 
the arch is a potential ridge of material. Hill et al.l (1201 ih define 
a 'ridge' as a high-column density (-lO^^cm"^) super-critical 
filament dominating its environment that is capable of form- 
ing high-mass stars. The most prominent part of the Northern 
Eilament has a column density above l.OxlO^^cm"^ over ~ 2pc. 
It is not as dominant, i.e., containing most of t he material in the 
vicinity, as the Vela C ridge or DR21 ridge (iHill et al.l l20 111 
Hennemann et al., in prep, respe ctively). This ridge and its s tar 
formation content is discussed by lWhite et al.l (12012. in prepI) . 

4. Impact of the NGC 661 1 cluster on the M1 6 cloud 

The NGC 661 1 cluster is responsible for most of the mechanical 
and radiative energy incident on the molecular cloud complex 
around the Eagle Nebula (Ml 6). Here we assess the dust heating 
eff'ects in view of constraining the geometry of the complex. 

4.1. The cold Eastern Filament 

The Eastern Eilament is a cold prominent filament running east 
of the NGC 661 1 cluster in the column density and dust temper- 
ature maps (Eigs.[2l[3]and|4]). This southern portion of this fila- 
ment houses many 70 yum sources (see Eig. lA.ll top-left), that we 
interpret as embedded protostars. In order to determine how the 
presence of the OB cluster in Ml 6, in particular its heating eff'ect, 
has impacted the Eastern Eilament, column density and dust tem- 
perature profiles were taken along the filament. These profiles 
were derived using the crest points identified by the DisPerSE 
algorithm (see Section 13.21) , and are presented in Eig. |5] (left) 
as a function of projected distance from the cluster. The dust 
temperature profile of the Eastern Eilament shows a decreasing 
temperature with increasing distance from the OB cluster. There 
are small fluctuations along the temperature profile, which are 
anti-correlated with column density (see the lower panel of Eig. 
[3 and thus density enhancements such as star-forming sources 
located within the filament. Around ~ 17pc from the OB clus- 
ter the temperature of the Eastern Eilament stabilises as it ap- 
proaches the southern edge of the map. At distances > 17 pc, the 
dust temperature in the Eastern Eilament increases with increas- 
ing distance as it approaches the Galactic Plane. The Eastern 



Table 1. The extent of the heating penetrating of NGC 661 1 on 
the two main filaments identified in the M16 complex, deter- 
mined from the dust temperature profiles (Eig. |5}. The param- 
eters of the column density and dust temperature derived from 
probability distribution functions (PDEs) for each sub-region in 
M16 (Section |B] and Eig. IB. II) as well as their heating distance 
and gradient are also given (Eig. ICIK The numbers given in 
parentheses in the last column is the value obtained when a de- 
screening has been applied (see section |4Jl) . The asterisk here 
(*) indicates that the heating gradient was taken over the steep- 
est portion (3-5 pc), see Eig. ICll 
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M16-NW 


yellow 


1.2 


17.2 
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Eilament has a ~ 3 K temperature gradient over its projected 
length of ~ 20pc. 

To remove the contributions of cold dense sources within the 
Eastern Eilament to the dust temperature we have a applied a 
de- screening to the data following the relatioi^^ by I Juvela et al.l 
(Eig. 1; 2011). This relation (green line on Eig. |5]) simulates 
the temperature structure for a filament with a homogeneous 
density crest, and and while it removes the eff'ects of den- 
sity enhancements above ref (where ref is taken here 
as 3.2 X 10^^ cm"^) including, cold dense fragments embedded 
within the filament, it does not remove the centrally heated pro- 
tostars (indicated as YSOs on Eig. [51). This column density filter- 
ing highlights the eff'ect of the OB cluster on the dust temperature 
of the Eastern Eilament as a whole, the overall trend of which is 
a decreasing dust temperature with increasing distance from the 
OB cluster up to a projected distance of ~ 17 pc. 

In contrast to the dust temperature, the mean column den- 
sity profile (Eig. [3 left) appears to remain rather constant with 
distance from the OB cluster, though it is subject to the afore- 
mentioned density enhancements. The strong (and wide) density 
fluctuation seen at ~ 20 pc reaching 7 x 10^^ cm"^, is located at 
approximately where the Eastern and Western Filaments appear 
to meet (see Section IT2l) . 

4.2. The Northern Filament 

The Northern Eilament originates at the eastern part of the 
NGC 661 1 cavity, run s aro und the cavity and across the poten- 
tial ridge (see Section l33l and White et ani2012. in prep) before 
continuing into the northern part of the map (see Fig. |4]l. The 
column density and dust temperature profiles of the Northern 
Eilament are given in Eig. [5] (right). The dust temperature of this 



^ Tcorrected = T,neasured + 4.4 log (N//2/ N//2 ref); whcrC N//2 ref IS the 

median column density of the filament. This relation, found in radiative 
transfer simulations of dense cores, is taken here as an ap proximation to 
that which should exist in Plummer-like (cf . Arzoumanian et al.ll2011h 
filaments. A paper is currently being prepa red on this topic which exam - 
ines such a relation in Herschel filaments (iDidelon et al.l2012. in prepI) . 
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Projected distance from the OB cluster [pc] Projected distance from tine OB cluster [pc] 

Fig. 5. The dust temperature (top) and column density (bottom) profiles along the Eastern (left) and Northern (right) Eilaments. The 
positions along the filaments are plotted as a function of their projected distance from the cluster. The NGC 661 1 cluster is located 
at '0' on the x-axis, centred at 18 18.8 -13 49. The green dashed line on the dust temperature (top) is t he dus t temperature after 
removing the eff'ects of cold sources within the filaments, following the relation of Juvela & Ysard (Eig. I: l201ll see footnote 6 and 
also Didelon et al., in prep, who examine this relation for Herschel filaments). The thick vertical lines seen on the profile of the 
Northern Eilament indicates structures at the same distance from the cluster, as the filament tends to curve around the cavity. 



filament remains rather constant (-20 K) around the NGC 661 1 
cluster, except for a temperature spike at ~ 3 pc corresponding to 
a position close to a well known Hii region (IRAS 18156-1343). 
At ~ 6 pc the dust temperature drops rapidly with increasing dis- 
tance from the NGC 661 1 cluster (5 K over 4pc), before stabil- 
ising at a near constant temperature of ~ 16 K. This temperature 
drop is greater than that seen in the Eastern Eilament (see Eig.[5l 
left). This diff'erence is (partly) because the column density pro- 
file is less homogeneous than the Eastern Eilament, with a large 
dip in the column density curve due to the presence of the ridge. 
While the dust temperature used here is a line of sight average 
for the filament, and the absolute values of the dust temperature 
are likely to be slightly overestimated as a result, the gradient 
identified in these Eastern and Northern filaments is still present 
and robust. 

Similarly to the Eastern Eilament, a de- screening was ap- 
plied to the Northern Eilament, with N//2 ref , the median col- 
umn density, taken as 4.4 x 10^^ cm"^. After applying this de- 
screening, the dust temperature of the Northern Eilament (green 
line Eig. [51 right) decreases steadily over a much greater dis- 
tance (5 K over 8pc). Even after de-screening, the drop in dust 
temperature is greater than that in the Eastern Eilament, which 
may result from the less homogeneous column density in the 
Northern Eilament or it might also suggest a larger relative in- 
clination along the line of sight, though spectral line (velocity) 
information is needed to support this idea. Note that at distances 
> 12pc the de-screening temperature is much lower (~ 2 K) than 
that without de- screening. This is likely due to the fact that there 
are very few dense cores in this portion of the filament (see Eig. 
O, and thus the de- screening does not completely apply here. 

4.3. NGC 6611 and its immediate surroundings 

In order to determine the extent of the heating from the OB clus- 
ter on low column density material (1 .5 x 10^^ cm"^ or lower, see 
Eig.O, PA-averaged radial profiles of both the dust temperature 
and the column density were taken covering all of the dense star- 



forming material in the region, as well as the cavity carved out 
by the NGC 6611 cluster where there is more diff'use material. 
These two radial profiles are presented in Eig. |6] with narrower 
radial profiles covering the interesting features in this region, 
e.g., the pillars, given in Eig. ICll The profiles, which take the 
median profile in the region rather than the average to limit the 
inclusion and bias of protostars, were centred at the position of 
the NGC 661 1 cluster. It is immediately clear (see Eig.O that the 
dust temperature is rather homogeneous close to the NGC 661 1 
cluster (~ 1.5- 2 pc) but then it drops with increasing distance 
from the cluster. The dust temperature profile covering most of 
the emission in the region (Segment 1, Eig. [6]) drops 2 K within 
~ 4 pc, and continues to decrease with increasing distance. The 
dust temperature profile of the cavity (Segment 2) indicates a 
homogeneous dust temperature close to the cluster (< 2pc), and 
then a steadily decreasing temperature with increasing distance 
from the cluster. The column density is anti-correlated with the 
dust temperature, a result which also applies to the smaller seg- 
ments presented in Eig. ICll Each of the segments presented in 
Eig. ICll display diff'erent dust temperature profiles, suggesting 
diff'erent column densities of material and/or three-dimensional 
eff'ects. 

The extent of heating by the NGC 6611 cluster on the sur- 
rounding environment, i.e., the penetration distance, is found to 
be 3 pc to 15 pc as determined from the dust temperature profiles 
(see the dust temperature gradient in Eigs. |5] [S] ICIK depend- 
ing on the region (see Table [TJ. It is not qualitatively diff'erent 
from the penetration measured for 3-6 x 10^^ cm"^ filaments 
(see S ection |4Jl and 14.21) . 



5. Discussion & Interpretation 

The Eagle Nebula (Ml 6) is arguably one of the finest exam- 
ples of high-mass star formation in the Galaxy, chiefly because 
of t he large number o f high-mass stars forming in this region 
(cf . iHillenbrand et al.l Il993i) and the remarkable visual efl'ects 
that this activity has had on the complex (e.g. the Pillars of 
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Fig. 6. Left: Column density map of NGC6611 with radial segments as indicated. Segment 1 (middle) contains all of the star 
formation activity in the region, and segment 2 (right) the cavity carved out by the NGC6611 cluster. These profiles have been 
averaged over azimuthal angle. The black solid line is the dust temperature profile, while the red dashed line is the column density 
profile. More detailed segments for each region inside segment 1 are given in Fig. lCll Cuts are the median profile in each segment 
to limit biases introduced by the presence of protostars. 



Creation, Spire, and the 'Eagle' in the optical). The Nebula's 
centre has largely been excavated by the radiation and power- 
ful stellar winds of the OB stars within the central NGC6611 
cluster. This radiation has shaped, and continues to shape, the 
molecular cloud material, possibly even triggering the formation 
of new stars and the dissipation of material. 

5.1. The (three-dimensional) structure of the complex 

The three-colour composite image of this region (Eig. [T) shows 
a very distinct temperature gradient running away from the cen- 
tre of the cavity carved out by the NGC 661 1 cluster. Those re- 
gions closest to the cluster are subject to heating by the OB stars, 
which is reflected in the higher dust temperatures (-1-3 -5 K) com- 
pared with those regions in the more remote (14-17 pc) parts of 
the cloud, where the heating no longer penetrates. The temper- 
ature gradient from the NGC 661 1 cluster is similar to that seen 
in Rosette ( Schneider et al. 2010) but over a smaller distance, 
with a clear drop in temperature within lOpc of the cluster in all 
directions (see Eig s. [5] and [6]) . 

The Ml 6 complex displays evidence of a large number of 
bright rims, which typically occur close to PDRs, especially in 
the vicinity of the NGC 661 1 cluster (see Eigs.O right and[3]and 
Section 13.21) Whilst a greater investigation into these potential 
sites of limb-brightening is warranted, we can use the informa- 
tion in these maps to infer the structure relative to the OB cluster 
(e.g., Minier et al., in prep). Those objects that exhibit signs of 
bright blue limb-brightening emission at their surface are likely 
directly receiving incident radiation from the NGC 661 1 cluster 
and are thus at the same distance as the cluster or behind it. There 
are clear pockets of hot dust (traced by 70 yum emission, shown 
in blue, in Eig.[T]) on the southern side of the NGC 661 1 cluster, 
around the pillars, and to the north of the cluster, as well as in- 
side the arch structure (also seen in Eig.|3j. This inner 10 pc shell 
is suggestive of a three-dimensional geometry where the bright 
70 yum objects, such as the Pillars of Creation, are illuminated by 
the NGC 6611 cluster and not extinguished by any cloud in its 
foreground. The location of the Spire and the ridge is less clear. 
On the other hand, objects such as the MYSO, which do not 
show evidence of bright rims, and thus incident radiation from 
the cluster, are clearly in front of the cluster, and closer to us 
along the line of sight. An interesting question is raised with re- 
spect to the origin of the arch structure (Section (33] and Eig. [4]). 



Is it actually one structure e.g., created by the expansion of radi- 
ation similar to that of a bubble, or is it a consequence of mor- 
phology? As noted in Section [33l two arches are in fact detected 
at this location and DisPerSE only detects the denser of these. 
Eigure [T] and Eig. [2] (right) suggest that the inner warmer arch 
is closer to the cluster, directly influenced by its heating, whilst 
the outer higher column density arch is likely above (i.e., closer 
to us along the line of sight) and away from the cluster. This 
suggests some sort of curved morphology, where the two arches 
are in fact two sides of a three dimensional shell-like or bubble- 
like structure, with the near (denser arch) and far (warmer arch) 
curved sides of the shell/bubble appearing as two arches in two 
dimensions. 

Two prominent filaments are detected in the Ml 6 complex, 
a cool dense Eastern Eilament (see Section |4Jl) which spans al- 
most 25 pc, and the Northern Eilament (see Section 14.21 ) which 
extends across the ridge, at the western side of the arch structure, 
and up to ~ 15 pc. Eigs.[2](left) and [3] indicate a clear branch to- 
wards the southern end of the Eastern Eilament, at ~ 15 pc from 
the origin (see Eig. [51 left), as discussed in Section |4Jl DisPerSE 
(Section 13.21) actually detects the much weaker (with respect to 
column density - Eig. [2]) western component of this branch (Eig. 
|4]), which extends up to the base of the Pillars of Creation. If 
we were to extrapolate or extend this western filament, it would 
continue to join the Northern Eilament at the approximate loca- 
tion of, and along the same line of sight as, the MYSO. In this 
scenario, the Eastern Eilament has remained largely undisturbed 
whilst the Western Eilament has undergone a dramatic transfor- 
mation. Once the home of the NGC 6611 cluster, it would have 
been destroyed as the cluster accumulated mass and formed the 
high-mass OB stars seen in existence today. The strong ionising 
winds from the OB stars themselves would have carved out a 
cavity in this filament, and eradicated any existence of the con- 
tinuation of this Western Eilament. Small grains associated with 
this Western Eilament are still present and emit strongly at 70 jim 
(Eig. lO top-left). 

5.2. The heating effect of the NGC 661 1 cluster 

The high-mass star cluster NGC 6611 has noticeably impacted, 
in terms of heating, the Ml 6 complex as a whole. The dust tem- 
perature profiles (Eigs. |5]and|6]) clearly show the effect and the 
extent of this impact, i.e., how much and how far it has af- 
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fected the temperature in all directions. Despite the fact that 
the dust temperature measured at the crest (Fig. |5} or aver- 
aged within segments (Fig. O obviously depends on the col- 
umn density of the material traced, the penetration distances do 
not seem very different from the diffuse emission of segment 1 
(< N//2 > 10^^ cm"^) to the high column density emission of the 
Northern Filament (<N//2 > 6xl0^^cm"^). Our analysis (see 
Figs. [5] and [6] and Table[T]) has shown that NGC6611 is able to 
heat the Ml 6 cloud to a projected distance of ~8 - 10 pc. The 
Ml 6 cloud displays -1.5 - 5 K temperature variations, corre- 
sponding to heating gradients of 0.25-0.6 K/pc, with an average 
value of -0.3 K/pc. The dust temperature profiles of the Eastern 
and Northern Filaments (Fig. O indicate that NGC 661 1 is able 
to heat, by 3-5 K, star-forming sites (> 10^ cm"^ dense cores) 
approximately up to 10 - 1 1 pc into each filament. While the dust 
temperature and column density gradients are robust, the abso- 
lute values of these quantities are subject to line of sight effects. 

It is clear from the radial profiles (Fig. |6) that the heating 
effect of NGC 661 1 is not the same in all directions. These dif- 
ferences would be expected if the NGC 661 1 cluster is in fact at 
the same distance as only some of these regions (as suggested 
in Section [5J1) . Alternatively, the ionisation and heating depths 
of the cluster is simply dependent on the amount of material in 
the way, such that high column density regions will block the 
radiation better than low column density regions in which the 
radiation more easily penetrates. 

Herschel observations have shown that interstellar filaments 
act as the birthplace of star formation. In Ml 6, we have shown 
that interstellar filaments can largely be influenced by heat- 
ing, such as that coming from the nearby high-mass NGC 6611 
cluster (e.g.. Fig. O. Such external heating is often not con- 
sidered when deriving the evolutionary status of a core from 
its far-infrared SEP . Evolu tion/age pro xies such as hsubmm/Lboi 
(lAndre et al.ll2000h or Ttoi (iMvers & L add 1993) for evolution- 
ary diagrams have been recently extrapolated to massive dense 
cores (e.g. Motte & Andre 2001) . These parameters should be 
taken with caution, since massive dense cores generally form 
in molecular cloud complexes associated with OB clusters. By 
ignoring this effect, T^oi will appear cooler at larger distances 
from the cluster, and hotter closer to the cluster. Similarly, the 
^submm/Lboi r^tio of a protostellar core will be larger further from 
the cluster than closer to it. Inferring evolutionary trends from 
these parameters, one would then naturally find more evolved 
protostellar dense cores closer to the cluster than further away 
from it, giving the (incorrect) impression of triggered star forma- 
tion. In order to avoid these biases, the effects of external heating 
on interstellar filaments should be considered when determining 
the evolutionary state of a star-forming core, especially in high- 
mass star-forming regions. 

As an example, a ~ 6 K dust temperature gradient has been 
found for massive dense cores forming in the Rosett e molecular 
cloud under the influence of the NGC 2244 cluster (iMotte et al.l 
Even after subtracting a possible 3 K gradient coming 
from the different temperatures of the cloud filaments inside 
which prestellar cores have formed and evolved into protostars, 
the remaining 3 K could still ar gue for an age t r end fo r these 
YSOs, as originally suggested by lSchneider et al.l (l2010h . 

Our Herschel data of Ml 6 do not provide any e videnc e 
for a supernova remnant, as suggested by Flagey et al] (l201ll) . 
Our maps (Fig. O right) indicate dust temperatures in Ml 6 
that ar e consistent with other high-mass star-formin g re- 
gions (Nguv Luong et al.ll201 lli fHifl et al. 2011). Flage v et al.l 
(I2OIII) . with their shorter wavelength Spitzer data, were more 
sensitive to hot dust which is being shattered by a supernova 



remnant. Based on this comparison alone, we are unable to sup- 
port or refute this claim based on our Herschel data, which are 
more sensitive to colder dust. Howeve r, young supernova rem- 
nant s such as SN 19 87 A dMatsuura et al.ll201 1) and Cassiopeia 
A (Barlo w et al .1 120 id) have been readily observed by Herschel, 
as have much older supernova remnants e.g. W44. SN 1987 A is 
seen as a point source with PACS and as extended emission at 
increasingly longer Herschel wavelengths, while Cas A is well 
traced at all Herschel wavelengths, as well as at 24yum with 
Spitzer. The fact that we don't see evidence of such a supernova 
remnant in our individual Herschel images, and that we see no 
evidence of impact on the column density distribution around the 
NGC 661 1 cluster, as well as the absence of emission seen with 
Spitzer fe.g. lFlagev et al.ll201 iHlndebetouw et al.ll2007l) does not 
favour the presence of a supernova remnant in this region. 
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Appendix A: Herschel images of M16 
Appendix B: Probability distribution functions 

To examine the dust temperature and column density structure in 
the entire Ml 6 region probability distribution functions (PDFs) 
were derived (Fig. IB. II ). Column density PDFs have been used 
by many authors to characterise cloud structure based on star for- 
mation activity ( Kainulainen et al. 2009; Federrath et al. 2010). 
The column density in Ml 6 spans almost two orders of mag- 
nitude (5.0 X 10^^ - 3.7x lO^^cm"^), whilst the dust tempera- 
ture ranges from 14.5 K to 23.5 K. The highest temperature in 
the maps is recorded directly west of the base of the Pillars of 
Creation and thus south of the NGC6611 cluster. In order to 
investigate the variation of PDF profiles within Ml 6, we segre- 
gated the complex into six sub-regions (see Fig.O left), accord- 
ing to the prominent structures or cloud features impacted by the 
NGC6611 cluster. These sub-regions include: the main compo- 
nent of material, corresponding to the Eagle's wings in optical 
images and the arch structure (hereafter referred to as the 'M16- 
N' sub-region; green in Fig. IB. IK the Pillars of Creation (blue). 
Spire (cyan), the MYSO detected by Indebetouw et al. (magenta 
l2007h . the eastern portion of the map (hereafter 'M16-E'; red), 
and the northern periphery (hereafter 'M16-NW'; yellow). The 
MYSO is treated as a separate sub-region due to its high column 
density tail, to prevent this MYSO from skewing the PDFs of 
other sub-regions. While the column density PDF of this MYSO 
should be treated with caution, its inclusion here allows compar- 
ison with the other sub-regions with respect to dust temperature. 



The column density PDFs are reliable for the larger re- 
gions (e.g., M16-N, M16-E, M16-NW) above Ay>10ma£l, 
whilst the smaller regions (e.g.. Pillars of Creation, Spire and 
the MYSO) are reliable only above Ay > 20 mag. We are only 
concerned with higher column density (>20mag) here. The 
M16-N, Spire and M16-E sub-regions have similar column den- 
sity slopes (around Ay > 30 mag). The M16-N and MYSO sub- 
region have high column density tails, similar to that seen in the 
Vela C Centre-Ridge sub-region (Hill et al. 20 ij), and consisten t 
with regulation by self-gravity (Balle steros-Paredes et al.|[2QTl]) . 

The Pillars of Creation sub-region has a column density 
PDF that resembles a log-normal shape which has been seen 
in clouds which are initially quiescen t and where turbulence i s 
likely the main shaping mechanism (iKainulainen et al.l l2009l) . 
Note, however that feedback eff'ects, such as those from a nearby 
OB cluster, have yet to be factored into turbulence simulations 
from which model PDFs are drawn. The suggest ion of turbu- 
lence here is consistent with the observations of Poundl (Il998l 
^^CO, ^^CO and C^^O) who found that the Pillars of Creation 
are not gravitationally bound and are being ripped apart. The 
sub-region hosting the Pillars of Creation (Fig. JB.ll blue) attains 
higher dust temperatures than the other sub-regions, including 
that of the Spitzer identified MYSO. This diff'erence shows that 
the NGC6611 cluster has had more of an impact on the cloud 
temperature than the current star formation process. The Ml 6-N 
and Spire sub-regions have similar dust temperature PDF slopes 
(> 21 K), which likely reflects their similar distance from the 
NGC 661 1 cluster, though the Ml 6-N sub-region contains colder 
gas in the northern part of the region, including the ridge itself. 
The M16-NW and M16-E sub-regions span only small ranges 
in temperature. The median column density and dust tempera- 
ture for each sub-region, determined from the PDFs, are given 
in Tableffl 

It is interesting to note that while the M16-E and Spire sub- 
regions have similar column density PDF profiles, they display 
quite diff'erent dust temperature PDF profiles. This suggests that 
they are both a part of the same cloud structure, but that the 
Spire is closer to and thus more readily externally heated by the 
NGC 661 1 cluster than M16-E. The three sub-regions closest to 
the OB cluster (Pillars of Creation, Spire, MYSO) have similar 
dust temperature profiles, but quite diff'erent column densities. 

There is no clear trend in the column density PDFs to sug- 
gest or represent an eff'ect of the winds or radiation from the 
NGC 6611 cluster on the sub-regions. The median temperature 
of each region correlates closely with the projected distance 
from the cluster, suggesting no strongly unbalanced projection 
eff'ects for the cloud surrounding the NGC 661 1 cluster. 

The probability distribution functions (Fig. IB. II ) confirm that 
the column density in each sub-region in Ml 6 is indeed dif- 
ferent, which influences the extent of the heating efl'ect of the 
NGC 6611 cluster on each sub-region, as seen in the dust tem- 
perature PDFs. At zero-order and above Ay > 15 mag there is no 
difl'erence seen in the structure of the cloud with distance from 
the NGC 661 1 cluster. That is, the column density PDFs of sub- 
regions near to the cluster, and thus heavily impacted, and those 
in more remote portions of the map, display similar column den- 
sity PDF slopes. 

It is also worth noting, that though normalised, the Pillars 
of Creation, Spire and MYSO sub-regions are comprised of a 
small number of pixels compared with the other sub-regions in 
the complex, which could thus afl'ect their interpretation. 



^ Where N//^ = 1 Ay x lO^i cm-^ mag-^ dBohlin et al.ll 19781) . 
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Fig. A.l. The 5 Herschel wavebands of Ml 6, from top 
to bottom, 70, 160, 250, 350 and 500 yum. Note also the 
extended PACS coverage compared with SPIRE (see 
Section [21). The observed level of cirrus noise in Ml 6 
is ~ 20mJy and 200 mJy dX A= 70 and 160yum, re- 
spectively and 2, 1 and 1 Jy at the SPIRE 250, 350 
and 500 yum bands. At shorter wavelengths only those 
warmer objects, such as protostars and Hii regions are 
seen. At longer wavelengths Herschel detects cold, 
deeply embedded filaments and the progenitors of high- 
mass stars. 
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Fig.B.l. PDFs of column density (top) and temperature (bot- 
tom) with the sub-regions as defined on Fig. [J] Only the high 
column density material that we are interested in > 10 mag is 
shown. 



Appendix C: Detailed radial segments around 
NGC6611 

The radial profiles of more detailed segments (than Fig. O 
around the NGC 6611 cluster, including the Pillars of Creation, 
can be found in Fig. ICll 
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Fig. CI. Upper left: Column density map of NGC6611 with radial segments as indicated. Upper-middle to Lower-right: Dust 
temperature (red) and column density (green) PA-averaged radial profiles in each segment. Segments are arranged clockwise starting 
with the Pillars of Creation (upper-middle) to the cavity region (lower-right) as identified by numbers. The profiles are the median 
values with projected distance in each segment to avoid biases by the presence of protostars. The radial lengths of each segment 
were defined by a drop in column density (as seen in the upper-left panel). 
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